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ABSTRACT 

The (V/Knax) of the cosmological X-ray flashes detected by Wide Field Cam- 
eras on BeppoSAX is calculated theoretically in a simple jet model. The total 
emission energy from the jet is assumed to be constant. We find that if the jet 
opening half-angle is smaller than 0.03 radian, the theoretical (V^/Knax) for fixed 
opening half-angle is less than ~ 0.4, which is consistent with the recently re- 
ported observational value of 0.27 ± 0.16 at the 1 a level. This suggests that the 
off-axis GRB jet with the small opening half-angle at the cosmological distance 
can be identified as the cosmological X-ray flash. 

Subject headings: gamma rays: bursts — gamma rays: theory 
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1. INTRODUCTION 

The X-ray flash (XRF) is a class of X-ray transients (Heise et al. 2001, see also Barraud 
et al. 2003). Some properties of XRFs, such as the observed event rate, the duration and the 
isotropic distribution, are similar to that of Gamma- Ray Bursts (GRBs), while the spectral 
hardness of XRFs characterized by the peak flux ratio, the fiuence ratio and the photon 
index is softer than that of GRBs. This class represents a large portion of the whole GRB 
population. Recently, the observational value of (V/Vmax), which is the measure of the 
homogeneity of spatial distribution (Schmidt, Higdon, & Hueter 1988; see also Chang & Yi 
2001; Kim, Chang, & Yi 2001), has been updated from 0.56 ±0.12 (J. Heise 2000, talk given 
in 2nd workshop Gamma- Ray Bursts in the Afterglow Era) to 0.27 ± 0.16 (J. Heise 2002, 
talk given in 3rd workshop Gamma- Ray Bursts in the Afterglow Era). The updated value 
of (V/Knax) suggests that XRFs take place at a cosmological distance. 

Various models accounting for the nature of the XRFs have been proposed (Yamazaki, 
loka, & Nakamura 2002; Heise et al. 2001; Dermer, Chiang, & Bottcher 1999; Huang, Dai, 
& Lu 2002; Mcszaros et al. 2002; Mochkovitch et al. 2003; Daigne & Mochkovitch 2003). 
Heise et al. (2001) proposed that XRFs could be GRBs at high redshift. The redshifts of 
XRF011030 and XRF020427 have an upper limit oi z < 3.5 (Bloom et al. 2003; J. Heise 
2002, talk given in 3rd workshop Gamma- Ray Bursts in the Afterglow Era), which suggests 
that XRFs take place at not so high redshift but the same as that of GRBs. The photosphere- 
dominated fireball model may account the nature of the XRFs with peak energy Ep more than 
~ 20 keV (Meszaros et al. 2002, Ramirez- Ruiz & Lloyd- Ronning 2002). However, further 
considerations are needed to explain the event with Ep ~ of approximately a few keV, such as 
XRF020427 (L. Amati 2002, talk given in 3rd workshop Gamma-Ray Bursts in the Afterglow 
Era), XRF020903 and XRFO 10213 (N. Kawai 2002, talk given in 3rd workshop Gamma-Ray 
Bursts in the Afterglow Era). The models with small Lorentz factors, such as the dirty 
fireball model (Dermer, Chang, & Bottcher 1999; Huang, Dai, & Lu 2002) or the structured- 
jet model (Rossi, Lazzati, & Rees 2001; Woosley et al. 2002; Zhang & Meszaros 2002a) also 
have possibilities to explain the properties of the XRF, with the implicit assumption that 
the XRF arises not from internal shocks (Zhang & Meszaros 2002b). The models for internal 
shocks with small contrast of high Lorents factors might be the origin of XRFs (Mochkovitch 
et al. 2003; Daigne & Mochkovitch 2003). 

For the other possibility, we have studied the off-axis jet model and proposed that if we 

observe the GRB jet with a large viewing angle, it looks like an XRF (Yamazaki, loka, & 
Nakamura 2002, Yamazaki, loka, & Nakamura 2003). In Yamazaki, loka, & Nakamura (2002) 
the value of the jet opening half-angle was adopted as A^^ = 0.1. In this model the distance 
to the farthest XRF ever detected is about 2 Gpc {z ~ 0.4) so that the cosmological effect 
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is small and (y/T4iax) ~ 0.5. Recent observations suggest that GRBs with relatively small 
opening angle exist, while the distribution of is not yet clear (Panaitescu & Kumar 2002). 

If we assume the total emission energy to be constant as in the previous paper, the intrinsic 
luminosity is larger for the smaller opening half-angle. Such GRBs at the cosmological 
distance observed from off-axis viewing angle may be seen as XRFs and {V/V-caax) is expected 
to be smaller than 0.5. 

In this paper, we will show that our off-axis model has a possibility of accounting for the 
observational value of (l^/Vmax) if we change some of the model parameters in the previous 
paper (Yamazaki, loka, & Nakamura 2002). This paper is organized as follows. In § 2, we 
describe a simple jet model including the effect of cosmological expansion. We assume the 
uniform jet with sharp edges. Although one may consider the structured jet motivated by the 
simulation of the coUapsar model, we cannot conclude, both observationally and theoretically, 
which model is preferable. The (y/Kiax) for the XRFs detected by the Wide Field Cameras 
(WFCs) on BeppoSAXis calculated in § 3. Section 4 is devoted to a discussion. Throughout 
this paper, we adopt the flat universe with = 0.3, = 0.7 and h — 0.65. 



2. EMISSION MODEL OF X-RAY FLASHES 

We consider a simple jet model of XRFs in the previous papers (Yamazaki, loka, & 
Nakamura 2002; loka & Nakamura 2001) taking into account the cosmological effect. A 

general formula to calculate the observed flux from an optically thin material is derived by 
Granot, Piran, & Sari (1999) and Woods & Loeb (1999). Here we adopt their formulations 
and notations. Let us use a spherical coordinate system r = (r, ^, 0) in the central engine 
frame, where the = axis points toward the detector and the central engine is located at 
the origin. Consider a photon emitted at time t and place r in the central engine frame. It 
will reach the detector at a time T given by 

T = (1 + z)T, ^{l + z){t- Tiijc) , (1) 

where \i = cos 9 and z is the cosmological redshift of the source; T = was chosen as the 
time of arrival at the detector of a photon emitted at the origin ai t — 0. Then the observed 
flux at the observed time T and observed frequency u, measured in erg s""^ cm~^ Hz"^, is 
given by 

jl,{n'^,r,T, + rfi/c) 



/?/i)2 

where d^, fl'^ and j^, are the luminosity distance to the source, the direction towards the 
detector measured in the frame comoving with the jet (comoving frame), and the comoving 
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frame emissivity in units of ergs s ^ cm ^ Hz ^ sr ^. The frequency v' , which is measured 
in the comoving frame, is given by 

v' = i/,7(l - = (1 + z)v-^{l - Pii) . (3) 

We adopt an instantaneous emission of infinitesimally thin shell aX, t — and r — tq. 
If the emission is isotropic in the comoving frame, the emissivity has a functional form of 



j'An'a,r,t) = Aof{v')5{t-t,)5{T-T,)H{M-\0~e,\) 



xH 



COS0 



cos A9 — COS 6y COS 6 
sin 9„ sin 9 



(4) 



where /(z^') and Aq represent the spectral shape and the amplitude, respectively. and 
6v are the jet opening half- angle and the viewing angle that the axis of the emission cone 
makes with the ^ = axis. The delta functions describe an instantaneous emission at t = to 
and r = tq, and H{x) is the Hcaviside step function which describes that the emission is 
inside the cone. Then the observed flux of a single pulse is given by 

^ ^ 2(1 + z)roc^'Ao M{T)f ^7(1 - /3cosg(T))] 

' dl mi-(3cose{T))f ' ^ ^ 

where, 1 - Pcos9{T) = (1 + z)-^{cP/ro){T - Tq) and Tq = (1 + z){to - ro/c/3). For > 



cos Ag-cos 0{T) cos t 



and < e{T) < AO - 9^, A0(T) = tt, otherwise A(f){T) = cos-^ Msine{T) 
< Ae, e{T) varies from to + A^ while from 9^ - Ad to + A^ for 9^ > Ad. 
latter case, A0(T) = for 9{T) — 9^ — AO. Pulse-starting and ending time are given by 



. For 
n the 



Tstart = To + (l + z)(ro/c/?)(l-/3cos(max[0,^^,- A^])), 
Tend = To + (l + ;2)(ro/c/?)(l-/3cos(^, + A^)). (6) 

The observed spectrum of GRBs is well approximated by the Band spectrum (Band et 
al. 1993). In order to have a spectral shape similar to the Band spectrum, we adopt the 
following form of the spectrum in the comoving frame. 



where as (Pb) is the low (high) energy power law index. Equations (5) and (7) are the basic 
equations to calculate the flux of a single pulse. 

In order to study the dependence on the viewing angle 9y and the jet opening half-angle 
A^, we fix the other parameters as ub — —1, Pb — —3, ^i'q — 200 keV, ro/c/37^ = 10 s, and 
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7 = 100 (Preece et al. 2000). We fix the amplitude Aq so that the isotropic 7-ray energy 
Eiso = 47rd|(l + z)-'^S{20 - 2000 keV) satisfies 

^^Ei,, = 5 X lO^^ergs, (8) 

when = and ^ = 1 (Frail ct al. 2001). Here 5(z/i - z/2) = /^'^"^^ F(T; z/i - z/2)rfT is 
the fluence in the energy range Ui — V2 and F{T] v\ — V2) = J^^ Fy{T)dv is the flux in the 
same energy range. The values of Aq for different opening angles are summarized in Table 1. 
When the jet opening half-angle becomes smaller, Aq becomes larger. 



3. CALCULATION OF {V/V^^) 

The formalism to calculate (V/Knax) is given by Mao & Paczyhski (1992) and Piran 
(1992). In our case the absolute luminosity and spectrum depend on 6^ and /S.6. The 
observed- integrated number count from the source at redshift z with the jet opening half- 
angle and the viewing angle 6^ is given by 

c(A0,e,..)=r"'iTri.^^!^4M^. (9) 

Let Zmin and Zmax be the minimum and the maximum redshift of the XRF for given A^ 

and 9„. In determining Zmin and Zmax-, we should note that the operational definition of 
the XRF detected by BeppoSAX is the fast X-ray transient with duration less than ~ 10^ 
seconds which is detected by WFCs and not detected by the Gamma-Ray Burst Monitor 
(GRBM). Therefore, if the sources are nearby such that z < Zmin, they are observed as 
GRBs because the observed fiuence in the 7-ray band (40 — 700 keV) becomes larger than 
the hmiting sensitivity of GRBM (~ 3 x 10~^ergs cm~^). If the sources are too far such that 
z > Zmax, they cannot be observed by WFCs, which have an observation band of 2 — 28 keV 
and a hmiting sensitivity of about 4 x 10~^ergs cm~^. Although the detection conditions of 
instruments vary with many factors of each event, such as the duration, the spectral index, 
or the peak photon energy (Band 2002), we adopt very simple criteria here. As shown in 
Figure 1, both Zmax and Zmin depend on 9^ and A^. 

For given z^ax, we can calculate the minimum integrated count Cmin — C{A9, 9^, Zmax)- 
Then, (V^/Kiax) for given A^ and 9^ is calculated as 

/_^\ j::::[ciz)/Cm.n]-'/^viz)dz 

WaxLe,e. i::::T^(')dz ' ^ ^ 
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where ^{z) is given as 



v(z)=^4j-^y ^(-^]. (11) 

^ ^ 1 + z \l + zj dz\l + z' ^ ^ 



The luminosity distance is given by 

dL{z) = f (1 + z) f J'' (12) 

The function n{z) is the comoving GRB rate density. We assume that n{z) is proportional 
to the comoving rate densities of the star formation nsriz) in the following form (Madau & 
Pozzetti 2000; Porciani & Madau 2001) 

nsF{z) = QAQh ^^P(^-^^) V^^ + Q^(l + ^)-3 Mo yr"^ Mpc-^ (13) 
exp(3.82;) + 45 

When we consider the case of z^nin < z^ax <^ 1, the cosmological effect can be neglected, 
so that "D oc and C/Cmin ~ {z/zmax)'"^- Then, Eq. (10) becomes 

^ 0.5 [1 + [Zmin/Zmaxf] . (14) 



Therefore, Zmin < Zmax < 1 imphes (l^/Kiax)A0,e„ ~ 0.5. 
Next, we integrate Eq. (10) over 9^ as 



Vmax/M JW{9,)de, 



(15) 



where W{6v) is the weight function which is the product of the solid angle factor and the 
volume factor: 

V{z) dz. (16) 

- min 

The results of the numerical integration are summarized in Table 1. For each A^, Zjnax{(^v,p) 
{zmin{(^v,p)) Hicans the maximum (minimum) redshift at which W{9v) takes the maximum 
value. If we take the jet opening half-angle as < 0.03, (T^/Kiax)A0 is smaller than ~ 0.4, 
which is consistent with the observational result at the 1 a level. 

Let us consider the behavior of (V^/Knax)Ae,6'„ in the case of the fixed opening half-angle 
A^ = 0.03. Recent analysis of the GRB afterglows shows that some GRBs have a jet with 
a opening half-angle of less than 0.05 radian, i.e., the smallest value of A^ is about 0.03 
(Panaitescu & Kumar, 2002). Let us decrease the viewing angle 9^ from a sufficiently large 
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value (7^^ ~ 20). As shown in Figure 1, both Zmax and Zmin increase monotonically, since 
the observed flux from the source increases due to the relativistic beaming effect. However, 
the behavior of (V^/V^ax)Ae,e„ is more comphcated since the function V^z) in equation (11) 
has the maximum value at z = Zp 1.5. We plot (V^/Knax)Ae,0„ as a function of in 
Figure 2. If 6^ is large enough Zmax is smaller than Zp. Then the cosmological effect is small, 
so that {V/Vma^)Ae,0y ~ 0.5. While if 9^ is small enough Zmin is larger than Zp. Then V^z) is 
a decreasing function in the range Zmin < z < z^ax , and the contribution of XRFs at smaller 
distance to (V^/Vniax)Ae,e^ is larger so that (V^/Knax)A6>,e^ is small. We note that the above 
argument does not depend on so much. 

The behavior of the weight function W{6y) is shown in Figure 3. When 6^ is large 
enough for Zmax to be smaller than Zp or when 9y is small enough, W{9y) is small since T^lz) 
or the solid angle factor are relatively small in the range Zmin < z < Zmax- One can see that 
in the case of = 0.03, W{9y) takes the maximum value at 7^^,p ~ 4.5, when Zmin ~ Zp. 

Above discussions show that the sources with the viewing angle ~ 0.05 at z ~ 1.5 
are the most frequent class of the XRFs in the population for the opening half-angle A^ = 
0.03. In other cases in which A^ is smaller than 0.03, the typical values are z ~ 1.5 and 
61^ ~ A6I + 0.02 since W{0y) takes maximum at B^^p ~ A6I + 0.02 (see Table 1). 

4. DISCUSSION 

We have calculated (V/Vmax) Ae for the emission from a simple jet model, and shown that 
when the jet opening half-angle A^ is smaller than about 0.03, (V^/Knax)Ae for the XRFs de- 
tected by WFCs on BeppoS AX is smaller than 0.4. The value of A^ ~ 0.03 has been obtained 
from the fitting of the afterglow fight curve (Panaitescu & Kumar 2002; Frail et al. 2001). 
Such a narrow jet which is inferred by afterglow observations is rare. However, following 
Eq. (1) of Frail et al. (2001), the jet break time is given by tj ~ 13(A^^/0.01)^'''^ minutes and 
so it requires fast localization to observe the jet break for a narrow jet. Therefore, at present, 
the small number of GRBs with small A^ may come from the observational selection effect. 
In the context of this scenario, we might be able to account for the fact that afterglows 
of XRFs have been rarely observed since the afterglow at a fixed time gets dimmer for an 
earlier break time. Furthermore, some "dark GRBs" might be such a small opening angle 
jet observed with an on-axis viewing angle for the same reason. 

We briefly comment on how the results obtained in this paper will depend on the Lorentz 
factor of the jet 7. We see that when 7 becomes large, {V/Vmax)Ae becomes small. For 
example, when we fix 7 = 200, we obtain (T^/Kiax)Ae=o.o3 — 0.39 and {V/Viaax)A9=o.i — 0.43. 



-8- 



This implies that the hmitation on can be relaxed. 

We can estimate the typical observed photon energy as huobs ~ (1 + z)~^6i'q, where 
— 7[1 — PcosO] and 9 — max{0, 9^ — A9} (Yamazaki, loka, & Nakamura 2002). Since 
7^1 and 9^, A9 <^ 1, we obtain 

hi^obs ~ (17) 

{1 + z)[l + i^9y] ^ ' 

In § 3, we have shown that for fixed < 0.03, the typical value of 9^ is ~ A^ + 0.02. 
Therefore, for the adopted parameters ^v'q — 200keV and the typical redshift z — Zp ~ 1.5, 
one can derive hugbs ~ 30 keV, which is the typical observed peak energy of the XRFs 
(Kippen et al. 2002). We can propose from our argument that the emissions from the jets 
with a small opening half-angle such as A^ < 0.03 are observed as XRFs when they are 
seen from off-axis viewing angle. If one can detect the afterglow of the XRF, which has 
the maximum flux at about several hours after the XRF, the fitting of light curve may give 
us the key information about the jet opening angle (Granot et al. 2002). Therefore, our 
theoretical model can be tested by the near-future observations. 

We can estimate the observed event rate of the XRF for fixed A^ as, 

RxRF,Ae^^J Wi9,)d9,. (18) 

In order to calculate -Rxrf, A6»j we consider the proportionality constant TZ = n{z)/nsF{z)- 
One can write approximately as 7^ = r^sM^k, where t^smq is the number of stars with masses 
M > 8Mq per unit mass. Since we assume all stars with masses M > 8Mq explode as core- 
collapse SNe, k represents the ratio of the number of XRF sources to that of core-collapse 
SNe. We adopt the value k — 1 x 10~^, which is derived by the result of Porciani & Madau 
(2001) combined with the effect of the solid angle factor (A^)^/2. Using a Salpeter initial 
mass function 0(M), we obtain r^sM^ = (/s^T® <^(^) ^^)/(/o^^^^° M0(M) dM) = 1.2 x 
IO-^Mq^ Then, in the case of A^ = 0.03, we derive Rxrf,A0 = 1 x 10^ events yr-^(7^/l x 
1O~^M0^), which is comparable to the observed event rate of the XRF (Heise et al. 2001). 
Note that the value of RxRF,Ae remains unchanged within a factor of 2 when we vary A^ 
from 0.01 to 0.07. 

When the jet opening half-angle A^ has a distribution /as, we integrate (V^/T4iax)A6» 
and i?xRF,A6> over the distribution of A^ as 

(y/V ) = / d{A9) fA0 RxRF, Ae (^/Kiax) Ae ^^^^ 
J d(A9) /ao RxRF,Ae 

^ _Jd{A9)fA9 RxRF, AO ,r,n\ 

"^^^^ - Jd{A9)fAe ' 
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respectively. We assume a power-low distribution as /a6i oc (A^) ^. When we adopt q — 4.54 
(Frail ct al. 2001), and integrate over A6' from 0.01 to 0.2 rad, we find (V^/Knax) = 0.36 
and -RxRF = 1 X 10^ events yr~^. These values mainly depend on the lower bound of the 
integration. For example, we obtain (V/Knax) = 0.43 and -Rxrf = 3 events yr~^ if the 
integration is done over from 0.03 to 0.2 rad (Note that we may let the value of -Rxrf 
be consistent with observed value by adjusting TV). Since the statistics of the observational 
data will increase in the near future owing to instruments such as HETE-2 and Swift, we 
will be able to say more than above discussion, including more accurate functional form of 
/a6> than that we have considered above, as well as the relation to the GRB event rate. 
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Fig. 1. — Maximum and minimum redshift of the XRF as a function of the viewing angle 
are shown in the case of = 0.03. The sohd hne (dashed hne) represents the maximum 
(minimum) redshift Zmax (zmin)- The jet emission is observed as the XRF if the source has 
a redshift z in the range z^in < z < z^ax (see text). The vertical dashed line represents 
Q^^lSB^ 0.03. 
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Fig. 2. — (V^/\4iax)Ae,6i„ for the XRF detected by the WFCs on BeppoSAX is shown as a 
function of the viewing angle jOy in the case of A9 = 0.03. 
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Y0V 

Fig. 3. — The weight function W{9y), which is the relative observed event rate, is shown as 
a function of the viewing angle 7^„. Note that the normalization of W{9v) is arbitrary. The 
vertical dashed line represents 9^ — A9 — 0.03. 
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Table 1. RESULTS OF CALCULATION FOR FIXED 







a b 








0.10 


0.84 


0.103 


2.8 


1.5 


0.46 


0.09 


1.0 


0.095 


2.9 


1.4 


0.45 


0.08 


1.3 


0.086 


3.0 


1.4 


0.44 


0.07 


1.7 


0.077 


3.1 


1.3 


0.44 


0.06 


2.3 


0.068 


3.3 


1.2 


0.44 


0.05 


3.4 


0.060 


3.5 


1.2 


0.44 


0.04 


5.2 


0.052 


3.6 


1.1 


0.43 


0.03 


9.3 


0.045 


3.8 


0.99 


0.40 


0.02 


22 


0.038 


4.0 


0.89 


0.38 


0.01 


109 


0.034 


4.1 


0.77 


0.35 



'^In units of erg cm ^ Hz ^. 

'^Thc viewing angle at which the weight function W{6y) 
takes the maximum value. 

^For the XRFs detected by WFCs on BeppoSAX. 



